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The phase transition to a deconfined phase is studied and the consequences in the formation of 
neutron stars are investigated. We use the quark-meson-coupling model for the hadron matter and 
the MIT Bag model for the quark matter in order to build the appropriate equations of state for 
'^f^ ' the hybrid stars. The properties of the stars are then calculated. The differences between unpaired 

' and color-flavor locked quark matter are discussed. 
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It is widely believed that hadronic matter undergoes a phase transition to quark matter at high 
^ ' densities and/or high temperatures. The high temperature limit is important in heavy ion collisions 

10 ! and/or cosmology, whereas the high baryon density behavior is important for the study of neutron, 

' hybrid and quark stars. 

\ After a gravitational colapse of a massive star takes place, a neutron star with practically zero 

temperature can be born few seconds after deleptonization. The correct calculations of the star 
properties depend on the appropriate equations of state (EOS) that describe its crust and interior. 
, The crust of the neutron star, where density is low is believed to be adequately described by hadronic 

matter. Its interior, however, where density is of the order of 5 ^ 10 times nuclear saturation density 
remains to be properly understood. Whether the central part of the star is composed of quark matter 
only, of mixed matter or of paired quark matter is one of the subjects of the present work. 
O \ In the present paper we are interested in building the EOS for mixed matter of quark and hadron 

^ ■ phases. We employ the quark-meson coupling model (QMC) 0, 0j including hyperons in order to 

describe the hadron phase. In the QMC model, baryons are described as a system of non-overlapping 
MIT bags which interact through the effective scalar and vector mean fields, very much in the same 
way as in the Walecka model (QHD) yL Many applications and extensions of the model have been 



^ , made in the last years |4L l5l I6ll7ll8l l9lllC 

■ While the QMC model shares many similarities with QHD-type models, it however offers new 

opportunities for studying nuclear matter properties. One of the most attractive aspects of the 
model is that different phases of hadronic matter, from very low to very high baryon densities and 
temperatures, can be described within the same underlying model, namely the MIT bag model. In 
the QMC, matter at low densities and temperatures is a system of nucleons interacting through 
meson fields, with quarks and gluons confined within MIT bag. For matter at very high density 
and/or temperature, one expects that baryons and mesons dissolve and the entire system of quarks 
and gluons becomes confined within a single, big MIT bag. Another important aspect of the QMC 
is that the internal structure of the nucleon is introduced explicitly. It is found that the EOS for 
infinite nuclear matter at zero temperature derived from the QMC model is much softer than the one 
obtained in the Walecka model Q]- Also, the QMC model nucleon effective mass lies in the range 
0.7 to 0.8 of the free nucleon mass, which agrees with results derived from non-relativistic analysis of 
scattering of neutrons from lead nuclei |llj| and is larger in comparison with Walecka model effective 
mass. Consequently, we expect that contrary to the NL3 and TMl parameterizations of the non linear 
Walecka model (NLWM) with hyperons, the mass of the nucleon does not become zero at densities 
p < lOpo, po being the density of nuclear matter at saturation. At finite temperature, there arises 
yet another difference between predictions of QMC and QHD, namely the behavior of the effective 
nucleon mass with the temperature at fixed baryon density. While in QHD-type models the nucleon 
mass always decreases with temperature, in the QMC it increases. The difference arises because of the 
explicit treatment of the internal structure of the nucleon in the QMC. When the bag is heated up, 
uark-antiquark pairs are excited in the interior of the bag, increasing the internal energy of the bag 
In what follows we consider the QMC with a constant bag constant, e. g. not dependent on the 
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scalar meson field as in the modified quark meson coupling models (MQMC) ^3 -Contrary to the 
QMC result, the bag radius increases with density for all MQMC models with a meson dependent bag 
constant [13 ■ For densities not much larger than nuclear matter saturation density po the bags start 
to overlap which implies a breakdown of the model. In the present work we are interested in describing 
properties of nuclear matter for densities which go beyond the saturation density and therefore have 
chosen to consider the QMC model. 

For the quark phase we have chosen to use both unpaired quark matter (UQM) described by the 
MIT bag model 12, 13. Ilal and paired quarks described by the color-flavor locked (CFL) phase. 
Recently many authors [la, [13j l20> UM liave discussed the possibility that the quark matter is 

in a color-superconducting phase, in which quarks near the Fermi surface are paired, forming Cooper 
pairs which condense and break the color gauge symmetry |22j . At sufficiently high density the favored 
phase is called CFL, in which quarks of all three colors and all three flavors are allowed to pair. 

Once the mixed EOS is built by enforcing appropriate Gibbs criteria and chemical equilibrium 
conditions, the properties of the stars are calculated and discussed. We have restricted ourselves to 
the investigation of neutron stars (T — 0) and hence, neutrino trapping is not an expected mechanism 
because it is important only during the first seconds after gravitational collapse of the core of a massive 
star [2l[23. 

In recent works |25l . [2^ , two of us have studied the properties of mixed stars whose equations of 
state were built with the non-linear Walecka model (NLWM) for the hadron matter [^ and the MIT 
Bag [T^ [TsL [I^ [Tsj and the Nambu-Jona-Lasinio [23 models for the unpaired quark matter. In the 
first work |25l | the effects of temperature were investigated and in the second one the consequences of 
imposing neutrino trapping at fixed entropies were studied. In both works, only in very special cases, 
the interior of the stars were made of quarks only. In general, a mixed phase of hadrons and quarks 
was favored. 

In what follows we compare the properties of neutron stars obtained within the QMC model with 
hyperons and with both quark models discussed above, namely, the onset of hyperons, mixed and 
quark phases, strangeness content, the maximum allowed masses and core composition. 

The present paper is organized as follows: in section II the QMC model with hyperons is reviewed. 
In sections III and IV the unpaired quark matter and the CFL phase are described. In section V the 
mixed phase is implemented and the results are shown and discussed in section VI. Finally, in the last 
section, the conclusions are drawn. 



II. THE QUARK-MESON COUPLING MODEL FOR HADRONIC MATTER 



In the QMC model, the nucleon in nuclear medium is assumed to be a static spherical MIT bag 
in which quarks interact with the scalar and vector fields, cr, w and p and these fields are treated as 
classical fields in the mean field approximation. The quark field, ipq{x), inside the bag then satisfies 
the equation of motion: 



1 



tpqix) = , q ^ u,d,s, 



(2.1) 



where is the current quark mass, and g^, and g'^ and denote the_ quark-meson coupling constants. 
The normalized ground state for a quark in the bag is given by P, 



■0g(r,t) ^ J\fgexp{-ieqt/RB) 
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where 

eq^nq + Kb {^gl uj -f ^g^r^pog 
with the normalization factor given by 



- Rb m* 



Rb m* 



= 2R%f,{xq) [nq{nq - 1) + i?Bm;/2] 



(2.2) 



(2.3) 



(2.4) 



where ^Iq = ^x'^ + (Rbiti*)'^, m* = m° — 5^ cr, Rb is the bag radius of the baryon, and Xq is the 

quark spinor. The quantities il>q, Cq, Pq, Nq, fig, m* all depend on the baryon considered. The bag 
eigenvalue, Xq, is determined by the boundary condition at the bag surface 



h{Xq) = (3qjl{Xq) 



(2.5) 
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The energy of a static bag describing baryon B consisting of three ground state quarks can be expressed 
as 

where Zb is a parameter which accounts for zero-point motion and Bb is the bag constant. The set 
of parameters used in the present work is given in tabled The effective mass of a nucleon bag at rest 
is taken to be 

M*b=eI''^. (2.7) 

The equiUbrium condition for the bag is obtained by minimizing the effective mass, Mg with respect 
to the bag radius 

-— ^ = . (2.8) 
For the QMC model, the equations of motion for the meson fields in uniform static matter are given 

by 



2Jb + 1 MUa) 



^0 - ^.9c.s (2./B + 1) 4/(67r') : (2-10) 



Y.3pBhB {2Jb + 1)4/(6^') • (2.11) 

B 

In the above equations Jb, I3B and kB are respectively the spin, isospin projection and the Fermi 
momentum of the baryon species B. The hyperon couplings are not relevant to the ground state prop- 
erties of nuclear matter, but information about them can be available from the levels in A hypernuclei 

M 

QaB = XaB 9c!N, QujB = X^B dujN, QpB = XpB QpN 

and x„B, Xi^B and XpB are equal to 1 for the nucleons and acquire different values in different pa- 
rameterizations for the other baryons. Note that the s-quark is unaffected by the sigma and omega 
mesons i.e. 5^ = gf, = . 
In Eq. ifrg)) we have 

g.BCBia) = = ^ Y: ^^^9lSB{a) (2.12) 

q—u,d 

where 

Sb{<j)^ dr^,V, = 7^775 uJj, J. ; (2.13) 

Jhag ^ ^qiS^q - 1) + RBm*/^ 

The total energy density and the pressure including the leptons can be obtained from the grand 
canonical potential and they read 

I22 I22 I22 
+ E '^-^^ ^"dk + A/^2(a)] + E ^ P ''"dk [k' + rnf] , (2.14) 

B ^ I ^ 

D I221I221I22 



2 

1 2Jb + 1 k* dk 1 1 /■'" k'^dk 



1/2 

'/J 



(2.15) 
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The lepton Fermi momenta are the positive real solutions of {k^ + iri^y^^ = fie and {k^ + m^y^^ = 
fJ-fi = Me- The equilibrium composition of the star is obtained by solving the set of Eqs. H2.9|) - 
(|2.1HI in conjunction with the charge neutrality condition (j2.16|l at a given total baryonic density 
p — X]b(2^-B + 1)^5/(6'''^); the baryon effective masses are obtained self-consistently in the bag 
model. 

For stars in which the strongly interacting particles are baryons, the composition is determined by 
the requirements of charge neutrality and /?-equilibrium conditions under the weak processes Bi — > 
B2 + I + Vi and B2 + I ^ Bi + vi. After deleptonization, the charge neutrality condition yields 

9tot = 5]te(2JB + l)4/(67r2) + qikf/{3n') = , (2.16) 

B l=e,fi 

where qs corresponds to the electric charge of baryon species B and qi corresponds to the electric 
charge of lepton species I. Since the time scale of a star is effectively infinite compared to the weak 
interaction time scale, weak interaction violates strangeness conservation. The strangeness quantum 
number is therefore not conserved in a star and the net strangeness is determined by the condition 
of /^-equilibrium which for baryon B is then given by /xb = 6b — qBfJ-e, where ^b is the chemical 
potential of baryon B and 5b its baryon number. Thus the chemical potential of any baryon can be 
obtained from the two independent chemical potentials /i„ and /ig of neutron and electron respectively. 

We start by fixing the free-space bag properties for the QMC model. For the bag radius Rn = 0.6, 
we first fixed the two unknowns Zj^ and Bf^ for nucleons. These are obtained by fitting the neuclon 
mass M = 939 MeV and enforcing the stability condition for the bag at free space. The values 
obtained are Zn = 3.98699 and B^J^ = 211.303 MeV for m„ = to^ = MeV and Zn = 4.00506 and 

1 /4 

= 210.854 MeV for m„ = nid = 5.5 MeV. We then fixed these bag values, Bb, for all baryons 
and the parameter Zb and i?B of the other baryons are obtained by reproducing their physical masses 
in free space and again enforcing the stability condition for their bags. The values obtained for Zb 
and Rb are displayed in table I for m„ = nid = MeV and for m„ = md = 5.5 MeV. Note that for a 
fixed bag value, the equilibrium condition in free space results in an increase of the bag radius and a 
decrease of the parameters Zb for the heavier baryons. 

Next we fit the quark- meson coupling constants g^, g^j — Sg^j and gp — for the nucleon to 
obtain the correct saturation properties of the nuclear matter, Eb = e/p — M = —15.7 MeV at 
p = pg = 0.15 fm-3, a,y„, = 32.5 MeV, K = 257 MeV and M* = 0.774Af. We have g$ = 5.957, 
gtjN = 8.981 and gpN — 8.651. We take the standard values for the meson masses, mo- — 550 MeV, 
= 783 MeV m'p = 770 MeV. 

For the meson-hyperon coupling constants we have opted for three sets discussed in the literature 
: set a) based on quark counting arguments we take x^b = x^ub = XpB ~ as in [S^; set b) 

according to [Sfl IsJ we choose the hyperon coupling constants constrained by the binding of the A 
hyperon in nuclear matter, hypernuclear levels and neutron star masses [xa = 0.7 and — Xp — 
0.783) and assume that the couplings to the S and S are equal to those of the A hyperon; set c) based 
on the 5f7(6) symmetry for the light quarks (u,(i) counting rule @ we take Xr,B — x^b — 2/3 and 
XpA = 0, Xps = 2,Xps = 1. 

In Fig. ^we have plotted the EOS obtained with the above parametrization of QMC and two choices 
of the hyperon couplings, set a) and b) together with the corresponding EOS obtained with NLWM 
with cubic and quartic sigma terms {K =300 MeV and M* = 0.7 M). Although the NLWM-EOS is 
harder at low densities, it becomes softer at higher energies after the onset of the hyperons. This fact 
has consequences on the behavior of the mixed phase and on the star properties. Moreover, one can 
see that different choices of the meson-hyperon parameters have a greater influence on the NLWM, 
for which the curves separate at lower densities and become more distant one from the other at high 
densities, than in the QMC. 



III. UNPAIRED QUARK MATTER EQUATION OF STATE 

The possible existence of quark matter in the core of neutron stars is an exciting possibility [T^ . 
Densities of these stars are expected to be high enough to force the hadron constituents or nucleons to 
overlap, thereby yielding quark matter. We take the quark matter equation of state as in Ref. [11 m 
in which u,d and s quark degrees of freedom are included in addition to electrons. Up and down 
quark masses are set to 5.5 MeV and the strange quark mass is taken to be 150 MeV. In chemical 
equilibrium pd = fJ-s = + Me ■ In terms of neutron and electric charge chemical potentials /i„ and 
/Lte, one has 

12 11 11 ,^ , 

M« = gMn - gMe, Md^gMn+gMe, Ps = ^^J'n + ^^J'e■ (3.1) 
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The pressure for quark flavor /, with f = u,d or s is 

= i E \^'fkM - 2-5™/) + l-^rnp 



where the Fermi momentum is fc/ = (/^^ — m^)^/^. 
For the leptons, the pressure reads 



^if + kj 



(3.2) 



Pi = 



—y 

37r2 ^ 



(3.3) 



The total pressure, including the bag constant B, which simulates confinement becomes 

P = Pt+Pg-B . (3.4) 

There are only two independent chemical potentials /Lt„ and /ig. /ig is adjusted so that the matter 
is electrically neutral, i.e. dP/dfi^ = 0. 



IV. COLOR-FLAVOR LOCKED QUARK PHASE 

In this section we study the equation of state taking into consideration a CFL quark paired phase. 
We treat the quark matter as a Fermi sea of free quarks with an additional contribution to the pressure 
arising from the formation of the CFL condensates. 

The CFL phase can be described with the help of the thermodynamical potential which reads 0| : 

^CFL{lJ-q,t^e) = riq„arfcs(Mg) + ^GB(A^9,Me) + ^i(Me), (4.1) 

where ^.q = /i„/3 and 



6 r 2 3 r 2 3AV^ 

^quarksi^^q) ^ — P dp{p ^ f^q) + — P dp^ p^ + - ^g) h B, 

Jo ■"■Jo 



(4.2) 



with m„ = nid set to zero, 



i. = 2A.,-V/^? + ^, (4.3) 



^Gsip-q, ^J■e) is a contribution from the Goldstone bosons arising due to the chiral symmetry breaking 
in the CFL phase [HEl: 



,2 N 2 



nGB{pq,^^e) = -lf^f^l (l-^j , (4.4) 
where 

(21-8 ln2)f,l 2 3A2 

•^'^ ^ 367r2 ' "^'^ ^ + '^'^)' y^-^) 

f2;(/ie) is the negative of expression (|3.3() . and the quark number densities are equal, i.e., 

Pu ^ Pd = Ps ^ 5 -■ (4.6) 

In the above expressions A, the gap parameter, is taken to be 100 MeV [l7| . 
The electric charge density carried by the pion condensate is given by 

QcFL = f^l^e(l-^y (4.7) 

In the above thermodynamic potential, we have neglected the contribution due to the kaon conden- 
sation which is an effect of order and thereby small compared to theA^^^ contribution to the 
thermodynamic potential for A ^ 100 MeV. 



6 



V. MIXED PHASE AND HYBRID STAR PROPERTIES 



We now consider the scenario of a mixed phase of hadronic and quark matter. In the mixed phase 
charge neutraUty is imposed globally i.e. the quark and hadronic phases are not neutral separately 
but rather, the system prefers to rearrange itself so that 

XPr + (l-x)pr + Pc = (5.1) 

where and p^^ are the charge density of quark and hadron phase, x is the volume fraction 
occupied by the quark phase, (1 — x) is the volume fraction occupied by the hadron phase and is 
the lepton charge density. As usual, the phase boundary of the coexistence region between the hadron 
and quark phase is determined by the Gibbs criteria. The critical pressure and critical neutron and 
electron chemical potentials are determined by the conditions. 



fiHP,i = P-QPA = fJ-i, i = n,e, Thp = Tqp, Php{p-hp,T) = PQp{fj.Qp,T), 

reflecting the needs of chemical, thermal and mechanical equilibrium, respectively. The energy density 
and the total baryon density in the mixed phase read: 

e = xe'^'' + {l-x)e"''+e\ (5.2) 



P = XP^^ + (1-X)P^^. (5.3) 

Notice that in all equations above the quark phase (QP) can be either the UQM or the CFL phase. 
The EOS for the mixed phase are then constructed. Once they are obtained, the properties of the 
neutron stars can be computed. The equations for the structure of a relativistic spherical and static 
star composed of a perfect fluid were derived from Einstein's equations by Oppenheimer and Volkoff 
They are 

dP _ G[s + P] [M + Airr^P] 

d7~~7 (r - 2GM) ' ^ ■ -* 



— — = 47rr^e, (5.5) 
dr 

with G as the gravitational constant and M{r) as the enclosed gravitational mass. We have used 
c = 1. Given an EOS, these equations can be integrated from the origin as an initial value problem for 
a given choice of the central energy density, (eo). The value of r {= R), where the pressure vanishes 
defines the surface of the star. We solve the above equations to study the structural properties of the 
star, using the EOS derived above. 



VI. RESULTS AND DISCUSSION 



In all figures shown, set a) for the meson-hyperon coupling constants was used, unless stated other- 
wise. We have omitted all results for the parameter sets b) and c) because they are very similar to the 
ones obtained with set a) . Actually, all differences appear only in the hadron phase at densities where 
the mixed phase is already the dominant one. It is worth emphasizing that, as stated in sections III 
and IV, the u and d quark masses are different in the UQM and CFL models. 

In figures[2t') and[2l3) the EOS obtained with the unpaired quark model and the color flavor locked 
phase are displayed for different values of the Bag pressure B. The onset of the mixed phase and quark 
pure phase occurs at lower denities for smaller values of B. This effect has already been discussed in 
|34j for the UQM description. A smaller value of B gives a softer EOS in the mixed phase because 
the onset of the mixed phase occurs at lower densities. However at higher densities, after the onset of 
the quark phase it becomes harder. The mixed phase shrinks with the decrease of the B parameter in 
both quark models. This fact can be also observed in tableni where the beginning and ending energy 
densities of the mixed phase are displayed in the last two columns. For the CFL, no mixed phase 
was found for B^/'*=180 MeV with iris = 150 MeV and A = 100 MeV, giving rise to a pure quark 
matter star. For the sake of comparison we have plotted the EOS for the same bag pressure for the 
UQM and CFL phase in figure 13 We include in the same figure the EOS obtained with NLWM plus 
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TABLE I: Bag constants for the baryons at the free space value B^^. The third and forth columns are 
obtained for B^'^ = 210.854 and the mass of the quarks taken as ruu = rud = 5.5 MeV and rua = 150 MeV. 
The last two columns are for Bg^'^ = 211.303 and the mass of the quarks taken as rrtu = irid ~ MeV and 
ms = 150 MeV 



Baryons 


Mb 


Zb 


Rb 


Zb 




Rb 


N 


939.0 


4.00506 


0.6 


3.98699 




0.6 


A 


1115.6 


3.69005 


0.62525 


3.68029 




0.62428 


E+ 


1189.3 


3.45577 


0.63977 


3.44628 




0.63870 


E" 


1192.5 


3.40386 


0.64038 


3.43600 




0.63931 


E" 


1197.4 


3.42970 


0.64038 


3.42024 




0.64024 




1314.9 


3.29260 


0.65336 


3.29188 




0.65182 




1321.3 


3.27173 


0.65455 


3.27105 




0.65301 


TABLE II: Mixed star properties 


model 


B^/^ (MeV) 


Mmax/MQ 


£0 (fm-*) 


£min (fm *) Emax (fm" 




QMC+UQM 


180 


1.41 


8.53 


1.26 


5.24 




QMC+UQM 


190 


1.58 


5.52 


1.63 


7.02 




QMC+UQM 


200 


1.73 


4.85 


2.05 


8.74 




QMC+UQM 


210.854 


1.85 


4.68 


2.73 


10.57 




QMC+CFL 


190 


1.32 


12.56 


1.35 


4.56 




QMC+CFL 


200 


1.49 


3.31 


1.92 


6.25 




QMC+CFL 


211.303 


1.76 


3.94 


2.66 


8.28 




NLWM+UQM 


180 


1.40 


7.38 


1.17 


4.62 




NLWM+UQM 


190 


1.64 


4.58 


1.81 


6.06 





UQM (short-dashed hne). One can see from this figure and tabled that the mixed phase in the CFL 
appears and ends at lower energy densities than in the UQM. This result has the same qualitatively 
behavior as the ones shown in Comparing QMC and NLWM with UQM one can see that they 
behave similarly at low densities. For e > 4.5fm-* the NLWM-EOS becomes softer in the mixed 
phase. This is due to the fact that at higher densities the NLWM with hyperons becomes softer than 
the QMC-EOS, Fig. d 

In figure 21 the particle population for the baryons, leptons and quarks are shown for the UQM and 
different bag pressures and for the CFL with one chosen Bag value. Hyperons only appear if the bag 
constant is very high, in our case B^/^ = 210.85 MeV. This value corresponds the Bb value of the bag 
used in the QMC model for the hadronic phase. A different behavior was obtained with NLWM in 
where the hyperons are present both for B^/^ = 180 and 190 MeV. As already discussed, both the 
mixed phases and the quark phase appear at lower densities for lower B values. For the CFL phase 
and _B^/* =200 MeV, the quarks appear at 2.67 nuclear saturation density. In for bag pressure 
of B^/^ =190 MeV, the quarks appear at 2.149 po. We can see that our results are compatible with the 
ones shown there, obtained within a NLWM formahsm with K = 240 MeV and M* = 0.78M. While 
in figures^]:), u, d and s quark populations come out different, in figure^Ji), because of the imposition 
of equal quark densities in the CFL model, they are forced to be equal. Also, the e~ population, 
which disappears at the onset of the quark phase in all the EOS studied, in the CFL model, goes to 
zero at a much lower density than in the UQM model, because the onset of a pure quark phase occurs 
at lower densities. 

As already discussed in [25I |2^ , the presence of strangeness in the core and crust of neutron and 
proto-neutron stars has important consequences in understanding some of their properties. In figure 
l^lwe show the strangeness fraction defined as 

rs=xrT ^i.^-X)r^'' (6.1) 

with 

OP ^ Pf. HP ^ Y.B \l?\PB 

6p 3p 

where qf is the strange charge of baryon B, for the different models discussed in this work. In all cases 
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the strangeness fraction rises steadily. If the UQM is used, at the onset of the pure quark phase it has 
reached 30 per cent of the baryonic matter. Although the amount of strangeness varies in the mixed 
phase, it is the same in the pure quark phase independently of the model used to describe the hadron 
phase. Nevertheless, if the CFL model is used for the quark phase, as a result of the equal quark 
densities imposition, the strangeness content reaches 1/3 of the total baryonic matter. Comparing 
QMC and NLWM with UQM, we conclude the strangeness content increases faster when the NLWM 
is used. This is due to the fact that in this model the hyperons also contribute. 

In tabled we show the values obtained for the maximum mass of a neutron star as function of the 
central density for some of the EOS studied in this work with UQM and CFL. Different bag parameters 
are used. We can see that the maximum mass of the star increases and its central energy decreases 
with increasing B. This result agrees with the fact that a larger B value corresponds to a harder 
EOS at high densities as shown in Fig. HI For B'^^^ = 180 MeV (QMC+UQM) and B^/'^ = 190 MeV 
(QMC+CFL), the central density of the star lies outside of the range of the mixed phase, an indication 
of a star with a quark core. In the last case we get a very high central density due to the very soft 
EOS this parametrization gives rise to. In fact the predicted maximum mass for a hybrid star within 
this parametrization is too low as discussed below. In all the other cases discussed the central density 
of the star is always within the mixed phase. We have also added some results obtained with the 
non-linear Walecka model instead of the QMC for the hadron phase for the sake of comparison. 
For the same B, the maximum mass is about the same, the central core is also made up of quarks 
and the mixed phase starts and finishes at lower densities. In |l7l| . the authors have found maximum 
masses around 1.6 Mq for maximum B values of B^^^ = 185 MeV. For these B values, our results 
come out at the same order. 

The radius of the maximum mass star is sensitive to the low density EOS. In order to calculate the 
radius and to plot it versus the star mass, we have used the results of Baym, Pethick and Sutherland 
for low baryonic densities. 

In Fig. Elthe mass of the family of stars obtained with QMC is plotted in terms of their radii for both 
quark models used and i?^/* ~ 190 and 200 MeV. We also include the family of stars obtained with 
NLWM plus UQM with S^/^ = 190 MeV. The radius of the stars within QMC is 11.95 Km (Si/* = 190 
MeV) and 12.47 Km (B^/^ = 2OO MeV) for UQM and 8.79 Km {B^/^ = 190 MeV) and 13.31 Km 
(51/4 = 200 MeV) for CFL. For the NLWM family of stars we get R = 12.53 Km (B^/-^ = 190 MeV). 
Some conclusions can be drawn. Comparing QMC and NLWM for B^/* ^ iqq jyjgy^ ^^le maximum 
mass of a stable star is similar, for both models. The fact that NLWM- Mq is larger shows that the 
main contribution to the star comes from the less dense regions, where the NLWM-EOS is harder than 
the QMC-EOS (see Fig. O. We consider now the families of stars obtained within QMC for both B 
values and using UQM and CFL. The quark contribution becomes more important for the smaller B 
values mainly for the CFL results. In particular, for B^/** = 190 MeV we get a quite small maximum 
Mq and a corresponding small radius. A maximum mass of 1.32 Mq is too small for accounting for 
the presently known radio-pulsar masses |3^ . even after corrections due to rotation. Similar numbers 
where obtained with the NLWM plus CFL in for slightly different values of B and m^. In fact, 
since these maximum mass stars have a small hadron exterior region the properties of the star is 
mainly determined by the quark model used. It is also clear from Fig. that the star properties 
are sensitive to the B value, in particular within the CFL formalism. Taking i?^/^ = 200 MeV the 
contribution of the mixed phase to stable stars is much smaller, restricted to the core of the star. In 
Fig. ^p) the properties of the stars which only contain hadronic matter do not coincide because we 
have taken ruu = md = 5.5 MeV for the QMC parametrization for the results QMC plus UQM and 
zero otherwise. 

A straightforward method of determining neutron stars properties is by measuring the gravitational 
redshift of spectral lines produced in neutron star photosphere which provides a direct constraint on 
the mass-to-radius ratio. Recently a redshift of 0.35 from three different transitions of the spectra of 
the X-ray binary EXO0748-676 was obtained in [37|. This redshift corresponds to MQ/R{km) = 0.15. 
In Fig. 6 we have added the line corresponding to this constraint and only the EOS for QMC plus 
CFL with B^/"'=190 MeV barely satisfies this constraint. In fact, the above constraint excludes all 
the EOS with hyperons, quarks or obtained within a relativistic mean-field approach, as compiled in 

M 

VII. CONCLUSIONS 

In the present paper we have studied the EOS for neutron stars using both the unpaired quark matter 
based on the MIT Bag model and the CFL phase for describing the quark phase and a relativistic 
mean- field quark- meson coupling description in which quarks interact via the exchange of a-, lo- and 
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p- mesons for the hadron phase. 

We have compared the results obtained within the QMC model with the ones obtained within the 
NLWM. For similar properties at nuclear saturation we conclude that contrary to NLWM, with the 
QMC the hypcrons only appear for a very high value of the bag constant. Also hyperons make the 
NLWM-EOS much softer than the QMC-EOS. For the hyperon meson couplings we have used three 
choices and verified that they did not have any effect on the onset of hyperons. The none appearance 
of hyperons affects the variation of the strangeness content of the star with density: except for the 
quark phase we get higher fractions with the NLWM. 

For the bag pressure parameter B we have used three different values, which produce different EOS 
and consequently the properties of the stars are dependent on them. Small B values give stable stars 
with a quark core. However, for a given gap constant A we obtain a phase transition to a deconfined 
CFL phase only if B is greater than a critical value. For lower values we get a EOS of pure quark 
matter, and pure quark matter stars. For A — 100 MeV we should have B^/^ > 185 MeV. For 
B^/-^ = 190 MeV and = 150 MeV, with QMC plus CFL we are not able to obtain stable stars 
with masses equal to most of radio-pulsar masses known. We have also shown that just one of the 
EOS we have studied satisfies the constraint imposed by the recently measured redshift of 0.35 from 
three different transitions of the spectra of the X-ray binary EXO0748-676 [37| . A more systematic 
study has to be done in order to verify whether a set of acceptable parameters for the QMC-I- MIT 
bag model gives a family of stars which contains the measured value. 

In 0, strange meson fields, namely the scalar meson field /o(975) and the vector meson field </)(1020), 
were also considered in order to reproduce the observed strongly attractive AA interaction. They have 
shown that the introduction of these strange mesons makes the EOS harder due to the repulsive effect 
of the 0(1020), meson. A harder EOS for the hadronic matter gives rise to a onset of the mixed phase 
at lower densities and a smaller mixed phase. The inclusion of these mesons and their influence on 
the properties of the stars are under investigation. 

It has been shown that the effect of temperature on the maximum mass of stable stars is small 
compared to the effect of neutrino trapping ^25„ ,26. 34] . Therefore, it would be interesting to include 
neutrino trapping even at T = MeV and impose Icptonic number conservation in the models used 
here and check the properties of the arising stars. 
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FIG. 1: EOS obtained with QMC and set a) (solid line), QMC and set b) (dot-dashed line), NLWM and set 
a) (dashed line) and NLWM and set b) (dotted line) including hyperons. 




FIG. 2: EOS obtained with QMC plus UQM (a) for B^/* 
MeV (long-dashed Une) and QMC plus CFL (b) for B^/*-- 
line) 



=180 (solid line), 200 (short-dashed line) and 210.85 
=200 (short-dashed Une) and 211 MeV (long-dashed 
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FIG. 3: EOS obtained with QMC plus UQM (solid line), QMC plus CFL (dashed line) and NLWM plus UQM 
(dotted line) with B^/*=190 MeV. Only here we have considered niu = rud = MeV in both QMC models to 
compare the equations of state. 
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(c) (d) 

FIG. 4: Particle fractions, Yi — pi/p, for i — baryons, leptons and quarks, obtained with the QMC+ UQM 
(a) for Bi/''=180 MeV, QMC+UQM (b) for B1/*=210.85 MeV QMC+UQM (c) for 6^/^=200 MeV, and 
QMC+CFL (d) for bV4=200 MeV 
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FIG. 5: Strangeness content obtained with QMC plus UQM (solid line), QMC plus CFL (dotted line) and 
NLWM plus UQM (dashed line) for B1/*=190 MeV. 
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FIG. 6: Neutron star mass versus radius obtained with the (a) QMC plus UQM (dotted line), QMC plus CFL 
(dash-dotted line) and NLWM plus UQM (long-dashed line) for 6^/^=190 MeV, (b) QMC plus UQM (dotted 
hne) and QMC plus CFL (dashed line) for Bi/''=200 MeV. The dots indicate the be ginning and end of the 
mixed phases. 



